Spin transport in interacting quantum wires and carbon nanotubes 
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We present a general formulation of spin-dependent transport through a clean one-dimensional inter- 
acting quantum wire or carbon nanotube, connected to non-collinear ferromagnets via tunnel junc- 
tions. We show that the low energy description of each junction is given by a conformally-invariant 
boundary condition representing exchange coupling, in addition to a pair of electron tunneling op- 
erators. The effects of the exchange coupling are strongly enhanced by interactions, leading to a 
dramatic suppression of spin accumulation. This is a direct signature of spin-charge separation in a 
Luttinger liquid. Furthermore, we demonstrate that magnetic polarization can lead to oscillations 
in the non-linear current-voltage relation. This phenomena is a surprising purely nonequilibrium 
effect due to backscattering interactions, which are thus dangerously marginally irrelevant in the 
repulsively-interacting Luttinger liquid. 



PACS: 71.10.Pm, 72.10.-d, 75. 70. Pa 

Recent studies on metal-ferromagnet hybrid systems 
have revealed new and interesting physics due to the in- 
terplay between the electronic charge and spin |J, e.g., 
the giant magnetoresistance effect Q . Following the ini- 
tial spin-injection proposal ||, the work of Johnson and 
Silsbee Q and subsequent advances have opened the way 
for the field of spintronics, where the electron spin is the 
central element for information storage and transport 
Spin-dependent transport plays an important role 
in quantum computation proposals, and has already led 
to new technological applications. In this context, a de- 
tailed understanding of transport through ferromagnetic- 
normal-ferromagnetic devices is both of fundamental and 
technological interest. In such structures, the current- 
voltage relation is predicted to sensitively depend on the 
relative angle 9 between the magnetization directions of 
the ferromagnets (FMs) 

Current theoretical models |(||7j are based on Fermi 
liquid theory, thereby ignoring the effect of interactions 
in the metal. As the inevitable miniaturization of spin- 
dependent devices proceeds, however, at least the inter- 
connects must ultimately reach the one-dimensional (ID) 
quantum limit, in which Fermi liquid theory breaks down 
||. This theoretically expected change from Fermi liq- 
uid to Luttinger liquid (LL) behavior drastically alters 
transport phenomena, as has recently been verified in 
experiments on charge conduction in carbon nanotubes 
H , which are nearly ideal ID quantum wires (QWs) [jlO| . 
Despite these developments, spin injection into a LL has 
received surprisingly little attention [jll). In this paper, 
we present a general low-energy theory for spin trans- 
port in a LL, which directly applies to nanotubes and 
semiconductor QWs [Q. We assume, as expected the- 
oretically [|l3| and recently observed experimentally fji|] 
for carbon nanotubes, that spin-orbit coupling in the LL 
is negligible. Its inclusion is, however, straightforward. 

Our analysis shows that spin transport in LLs is quai- 



itatively different both from charge transport in LLs and 
from Fermi liquid spin transport. We focus for concrete 
results on the case of an end-contacted quantum wire, 
and assume that the distance L between contacts is suf- 
ficiently long, max(V,T) 3> v/L, to ensure an incoherent 
stepwise transport mechanism through the tunnel barri- 
ers between each FM and the QW. (Here v is the Fermi 
velocity, and we put e = fog = h = 1.) A complete ana- 
lytic solution of this problem is contained in Eqs. (g) and 
( |r2| ) to (14). In contrast to charge transport, we find that 
spin conduction occurs not only through electron transfer 
but also exchange. This exchange effectively gives rise to 
a modification of the boundary conditions at the end of 
the LL, e.g., for the left contact, 



j r = n{Q)j L + J t 



unncl j 



(1) 



where Jl/r is the left/right moving spin current into/out 

of the contact, and Jtunnoi represents the effect of elec- 
tron transfer, see Eq. (|T2|). The effect of exchange 
coupling is given by the one-parameter SO(3) matrix 
TZ(Q) = exp(9r), where = ^a™^ 6 ^' an d r ™ is 
a unit vector in the direction of magnetization of the 
FM. Physically, represents the angle an incident spin 
in the LL precesses due to exchange interaction with the 
FM. Due to spin-charge separation in the LL, the ex- 
change contribution is not suppressed by the orthogo- 
nality catastrophe affecting the tunneling current, and 
therefore dominates the physics in many situations. This 
enhancement of the exchange current does not occur in a 
Fermi liquid, and its observation would provide a direct 
experimental signature of electron fractionalization. In 
addition to the novel physics arising at the contact, we 
find that a long ballistic QW exhibits a bulk precession 
of the magnetization, 



vd x M + d t J = bMx J, 



(2) 



1 



where M = Jr + Jl and J = Jr — Jl are the local mag- 
netization and current in the QW, respectively. Eq. (|^) 
leads to oscillations in the nonlinear current-voltage rela- 
tion. Remarkably, the latter is a purely non-equilibrium 
effect that arises from a marginally irrelevant backscat- 
tering interaction in the LL. The detailed character of 
these oscillations is also influenced by interactions. 

We now turn to the derivation of these results. In 
the incoherent limit, we may consider each contact sep- 
arately, as an initial system composed of two decoupled 
pieces, H — Hfu + Hqw- The FM {x < 0), described 
by Hpnj, is polarized along direction m, while the SU(2) 
invariant QW (x > 0) is described by Hq W . The SU(2) 
invariance guarantees the existence of a continuity equa- 
tion for spin density and current. At time t — > — oo, each 
half is assumed at equilibrium at its own chemical poten- 
tial, hfm and Hqw, and with a spin chemical potential 
h (see below) in the QW. We are interested in the steady 
state achieved at t = 0, long after the tunneling perturba- 
tion has been adiabatically turned on, H(t) = Ho + e St H' 
(5 — » + ). The calculation is non-trivial primarily due to 
its non-equilibrium nature: the system evolves according 
to H(t) while the initial states are distributed according 
to exp(— (3H ). Consider 

H' = F^W^ + ^W Jf F , (3) 

where F and \& are spin-1/2 Fermion annihilation opera- 
tors at x = for the FM and the QW, respectively. Em- 
ploying the projection operators u s = (l±m-<x)/2, the 2x 
2 tunneling matrix W takes the form W = J2 S t s u s , with 
spin-dependent transmission coefficients t s for spin quan- 
tization axis parallel to rh. The junction is then charac- 
terized by the conductance G — Gj + G± and the polar- 
ization P = (G T -G x )/G, where G u = {e 2 /h)\t tl \ 2 are 
the spin conductances J?]]. 

From Eq. (j^) and the spin continuity equation, the 
tunneling spin current across the junction is J = 
-i (ftffff* - &aW^F) j2. By defining H = H + 

HfmNfm + HqwNqw + h ■ Sqw, the standard pertur- 
bativc result can be rewritten as 



(J) 



Re ]T 

Q/37A 



dte 5t (Wa\ 



a(3 



x 



[^(0)^(0), mt)F X (t) 



(4) 



Thereby an intrinsically non-equilibrium expectation 
value is expressed in terms of an equilibrium aver- 
age using the shifted Hamiltonian Hq, where the non- 
equilibrium nature of the problem is fully encoded in 
the time dependent unitary matrix U(t) = exp[i(V + 
h ■ a/2)t], with V — [iQW — Hfm- A formula similar to 
Eq. (|j) can easily be written down for the charge cur- 
rent, I = i{F^W$> - ^W^F). Thus both charge and 
spin current can be calculated using equilibrium correla- 
tion functions. 



To proceed, we specify the Hamiltonians Hfm and 
Hqw- For energies well below the electronic bandwidth 
D, the F and 'J equilibrium correlators are identical for 
Hq and Ho j and moreover a non-interacting Fermi liquid 
model with constant density of states (DOS) applies to 
the leads. Because the lead couples to the QW only at 
x = 0, the difference in DOS for majority and minority 
spin carriers can be absorbed in a spatial rescaling of the 
Fermi fields of the FM and a suitable redefinition of the 
transmission coefficients (t s ) M. Then 



H 



FM 



dxfr z 8 x f , 



(5) 



where the spinor / = fop is indexed by b — (R, L) , de- 
scribing right- and left-moving modes of the FM, and 
by = (1\ 4) f° r the spin, with the boundary condition 
/fl(0) = /l(0). Here the Pauli matrix r 2 acts in the R/L 
space. Putting F = /(0), we see that F has SU{2) in- 
variant correlation functions. The low-energy description 
of the QW is an interacting LL model, 



QW 



= j dx {~i^vT z d x <ip + u(ipH) 2 } , (6) 



where i/>_r(0) = ^l(0) and "J" = ip(0). Only the forward- 
scattering interaction u is kept in Eq. (^). Alternatively, 
the exponent a > for tunneling into the end of the LL 
(x — 0) serves to measure the interaction strength ||. 
The SU (2) invariance of Eqs. (||) and (^J) implies 

Ft(O)*0(O),*t(i)F A (t)l) 9(-t) = 8 aX 5 M iC{-t) , 

where C(t) is the retarded Green's function of the oper- 
ator Ft^i (the choice of spin components is arbitrary). 
From Eq. (^) and the corresponding expression for /, it 
is then straightforward to obtain 



</> = -f£ 



(Pm + sh) Im C(V + hs/2 + iS) 



-Psm xhRe C(V + hs/2 + iS) 



(7) 



(7) = -G^(l + Psm ■ h) ImC(V + hs/2 + iS) , (8) 



where C(V) = Jdt C(t)e lVt . The terms involving 



Jq.(V, T) = GImC(V + iS) 
= GT(T/D) a smh(V/2T) 



(9) 



a V 

r 1 + - + i — 

' 2 2wT 



have a simple interpretation in terms of tunneling via 
Fermi's golden rule, as can be seen from the spectral 
representation of G. However, the appearance of ReG 
in Eq. M) indicates the presence of a physical process 



2 



other than tunneling. It can be shown |6|,|l3| that it cor- 
responds to a virtual process in which an electron near 
the Fermi energy in the QW hops into a state of the FM 
(which could be far from the Fermi energy) and hops 
back, thereby generating an exchange coupling. We thus 
include it from the start, H — » Hq + H' + H" , with 

H" = -Kfh-^a^/2 . (10) 

Since the FM possesses a non- vanishing average magne- 
tization, the spin operator in the FM may be replaced by 
this average to leading approximation. 

It is helpful to view both H' and H" in a renormal- 
ization group (RG) framework, as perturbations to a 
decoupled fixed point described by Hq. Standard ar- 
guments give the scaling dimension of both t s and K, 
At s = 1 + a/2 and Ak = 1. The scaling dimension Ak 
is not renormalizcd due to spin-charge separation in the 
QW. A simple calculation gives the RG scaling equations 

de\t a \ 2 (e) = -a\t B \ 2 , d e K(£)=c(\t r \ 2 -\t l \ 2 ) , 

where £ = hx(D/E), and c is a non-universal constant. 
Following the RG flow from the ultraviolet cutoff D down 
to energy E « max (T, V) <C D, we find \t s \ 2 (E) = 
\t s \ 2 (E/D) a , and 

K{E) =K + a- x cGP[l - (E/D) a ] » \t s \ 2 {E) . 

Therefore the tunneling spin current is much smaller than 
the exchange contribution. Neglecting the tunneling con- 
tribution completely, one still obtains a T-independent 
exchange spin current as T — » 0. 

This fact can be understood from a simple analogy to 
the Andreev current through a ballistic superconductor- 
normal-superconductor (SNS) junction p5[ . Let us con- 
sider a LL connected to two insulating FMs at x = 
and x = L, with rh x ml ^ 0. This is an equilibrium 
situation, which can be modeled using Eq. (^j) for the 
LL and two copies of Eq. (|h]) for the contacts to the 
FMs. The exchange interaction operates entirely within 
the spin sector of the LL due to spin-charge separation. 
Since the charge sector is decoupled, we are free to con- 
sider it at the non-interacting point, u = 0. Then the 
resulting fictitious charge boson and the physical spin 
boson can be combined and refermionized into a spin-ful 
Dirac fermion rj. Choosing arbitrary quantization axes 
rh = x and ml — cos{0)x + sin(0)y, it is instructive to 
perform the particle- hole transformation r/i —^rj\. This 
yields the Hamiltonian 

H„ = -iv [ dx rfr z d x r) - Re [KA(0)+K'A(L)e ie ] , 
Jo 

(11) 

where A(x) — r)^{x)rii(x). Equation ( |Tl| ) describes a 
ballistic SNS junction, and for phase twist < 9 < 2ir 



between the superconductors, supports an equilibrium 
current due to Andreev reflection. Since the Andreev 
current is vrfr z r], the original FM-LL-FM device indeed 
has a non-zero spin current J z . The analogy to a SNS 
junction also demonstrates that this current does not rely 
upon the incoherence of the two contacts. 

A more general perspective on the exchange coupling 
can be gained by viewing the low-energy physics entirely 
in terms of boundary operators and boundary conditions 
p6[ . For that purpose, we may make an arbitrary choice 
of short-scale physics, and let the exchange coupling act 
on right-movers slightly away from the junction. In this 
case, using the boundary condition ipi,(0) = ^r(O), the 
equations of motion for the spin currents can be inte- 
grated over the junction region to give the steady-state 
relation Jr(0 + ) = 1Z(®)Jl(0 + ) (the brackets denoting 
expectation values are omitted henceforth). The param- 
eter w K/v ultimately defines the "exchange coupling 
constant" of the low-energy theory. In principle, since 
the boundary exchange operator is exactly marginal, 
need not be small. Then the "bulk" spin currents are 
Jl = Jl(0 + ) and J R = Jr(0+) + J tU nnoi, and we obtain 
Eq. (Q), with the tunneling spin current 

/tunnel - ~ ^(^m + sh)l a {V + hs/2,T) . (12) 
s 

The term proportional to ReC* in Eq. (0) has been 
dropped, as its physical effects are included via the SO(3) 
rotation TZ. Since the magnetization far from the contact 
is M — \h with the LL spin susceptibility Xj one can 
then compute the spin current J for arbitrary exchange 
coupling 0. We arrive at the general result Jl7[ 

J = Sxh+ (1 - SVtunnel , (13) 

where S = (1Z—1) / (7^+1) is a real antisymmetric matrix. 
Similarly, the charge current is 

I = -^(1 + sm ■ h)l a (V + hs/2,T) . (14) 

s 

From Eqs. ( |l3| ) and (|lj), by exploiting spin and charge 
current conservation in order to obtain iiqw and h, one 
can then compute all transport properties in a given cir- 
cuit for arbitrary parameters |13| . 

We now specialize to a FM-LL-FM device with iden- 
tical contacts at T = and applied voltage V within 
v/L <C V <C D. For algebraic simplicity, we require 
P 2 <C (1 + a) 2 . For a tunneling contact, one expects 
< 1 (|, and Eq. (|||) then yields 

J= -(8x/2) rhxh + J tunnc i . 

Under these conditions, it is straightforward to obtain 
the ^-dependent FM-LL-FM current-voltage relation, 



3 



1(0) 



GV 



(V/D) a 1 - P 



tan 2 (0/2) 



tan 2 (6»/2) + r Q (V r ) 



(15) 



where Y a (V) = 1 + (6x/2) 2 (l + a)~ 2 (V/ D)- 2a . For 
a — > 0, this reproduces the result of Ref. @. Notably, 
unless the magnetizations of the FMs are anti-parallel 
(6 = 7r) or the exchange coupling vanishes (9 = 0), the 
spin accumulation effect Q|, in which the current is re- 
duced due to pile-up of spin in the QW, is strongly sup- 
pressed by the voltage dependence of Y a . Physically, this 
suppression is caused by the exchange coupling which is 
much more efficient in relaxing the injected spin polar- 
ization compared to the tunneling current. 

Finally, we turn to backscattering electron-electron in- 
teractions of the form 



H b = -b dx J L ■ Jr 



(16) 



For a carbon nanotube, with the lattice spacing a and the 
tube radius R, one may estimate b w ae 2 /R |10(| . Since 
Hb is marginally irrelevant in a single-channel QW, it is 
usually neglected in the LL model (||). Nevertheless, as is 
shown here, dynamical effects caused by ( filf ) can be im- 
portant. The equations of motion away from the contacts 
give Eq. (H) and vd x J + dtM — 0. In the steady state, 
we have conserved spin current J, and a bulk precession 
equation for M. Since M — the vector h must then 
precess around J. 

To get sizeable consequences, detailed analysis |l|] 
shows that it is essential to have small exchange cou- 
plings. For simplicity, we put = below. For the 
FM-LL-FM device, 



h(0) ■ h(L) = cos(A^) 



Aip = bJL/v 



(17) 



Since Aip oc L, precession will then always be significant 
for a sufficiently long QW. The computation of h(0) and 
h(L) leads to the self-consistency equation 



(1 - x 2 ) cos 



V 



cos(0/2) AV 



[V/DT = sin 2 (0/2) , (18) 



where AV = 8irv/[GPbLcos(9/2)]. Here solutions x = 
x n (n > 0) in the interval < x < cos(6>/2) have to be 
found. The current through the device is then 



In(V) 



GV 



(V/DY{l-P 2 [l-x 2 n (V)]). (19) 



The solution /„ corresponding to n full preces- 
sion periods exists only for voltages V > V n = 
J D[AT//D] 1 /( 1+Q )n 1 /( 1+Q ). Using typical parameters ap- 
propriate for a 1 (im long single-wall nanotube, one finds 
V\ rs 0.1 to 1 V. In general, the current- voltage rela- 
tion could then be multi-valued, where in the regime 
V n < V < V n +i, the solution I n (V) is expected to be re- 
alized. Hence the current-voltage relation becomes oscil- 
latory, with sawtooth-like oscillations. The observation 



of several periods could provide a direct and accurate 
measurement of a via the V n . 

To conclude, we have presented a general formalism 
for spin-dependent transport through interacting ID con- 
ductors. An experimental check of the theory should 
be possible by measuring the current- voltage relation for 
a ferromagnet-nanotube-ferromagnet device. The ap- 
proach is general enough to apply to numerous other 
problems. Several interesting extensions currently under 
investigation are the description of bulk- contacted wires, 
inclusion of the subband degree of freedom of single-wall 
nanotubes, LL to LL contacts, and the ballistic-diffusive 
crossover potentially relevant for multi-wall nanotubes. 
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